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1 INTRODUCTION 

• I— I . 

, The relationship between circum-nuclear star formation and 
5_j non-thermal nuclear activity is as yet poorly understood. 
Ci While the possibility that non-thermal activity may be en- 
t irely ascribed to massive stars and to supernova events 
( Terlevich fc Melnick 198^ ; ?) is challenged by several lines 
of evidence, it is unlikely that the onset of nuclear activity 
and strong nuclear or circum-nuclear star formation can be 
fully unrelated phenomena. Supernova may not be respon- 
sible for the optical spectrum and emission line of Seyfert 
1 galaxies, but reprocessed gas of supernova ejections may 
eventually be accreted by the central black hole, influenc- 
ing the opacity and hence the radiating properties of the 
accretion disk which is reputed to be one of the ultimate 
components of the AGN central engine. 

Core-collapse supernova from massive progenitors, al- 
beit rare events, are diagnostics of recent star formation. 
The aim of this paper is to estimate the ratio between 
the total number of Ib-Ic supernovae type II supernovae 
[5=N(Ib/c)/N(II)] discovered in a suitable sample of non- 
active (hereafter normal) galaxies, as well as to compare it 
to that of Seyfert galaxies. As we will show through this pa- 
per, the ratio S reflects metallicity, age, fraction of binary 
systems, and initial mass function (IMF) shape effects which 



ABSTRACT 

This paper estimates the relative frequency of different types of core-collapse super- 
nova, in terms of the ratio S between the number of type Ib-Ic and of type II super- 
nova. We estimate S independently for all normal and Seyfert galaxies whose radial 
velocity is <14000 km s~^, and which had at least one supernova event recorded in 
the Asiago catalogue from January 1986 to August 2000. We find that the ratio S is 
« 0.23±0.05 in normal galaxies. This value is consistent with constant star formation 
rate and with a Salpeter Initial Mass Function and average binary rate « 50%. On the 
contrary, Seyfert galaxies exceed the ratio S in normal galaxies by a factor w 4 at a 
confidence level ^ 2a. A caveat is that the numbers for Seyferts are still small (6 type 
Ib-Ic and 6 type II supernovae discovered as yet). Assumed real, this excess of type 
lb and Ic with respect to type II supernovae, may indicate a burst of star formation 
of young age (r ^ 20 Myr) , a high incidence of binary systems in the inner regions (r 
^ O.4R25) of Seyfert galaxies, or a top-loaded mass function. 
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are probably much different in normal and (at least in some) 
Seyfert galaxies. 



2 SAMPLE SELECTION & CRITERION 

We cross correlated the latest available Asiago Supernova 
Catalogue (last sup ernova recorded: 2000di on Aug. 23) 
(Barbon et al. 1993) with the 9**" edition of A Catalogue 
of Quasars and Acti ve Galactic Nuclei by M. - P. Veron- 
Cetty and P. Veron ([Veron-Cetty fc Veron200C|). All AGN 



hosting a supernova have Vr ^, 14, 000km s~ . A temporal 
restriction is needed sin ce type lb an d Ic supernova, flrst 
observed in the 1960s (Bertola 1964) were labeled as pe- 
culiar "type I" supernovae for more than two decades. The 
realization that they are a diffe rent class of objects came 



after two decades, around 1985 (Elias et al. 1985; Panagia, 



Sramek, fc Weiler 1986| ; porter fc Filippenko 198^ . We limit 



our counts to supernova discovered after Dec. 31, 1985. Only 
few identification of type Ib/c supernova discovered earlier 
than Jan. 1, 1986 are available (as a product of subsequent 
re- identification of older type I spectra). 

We define a sample of "normal" galaxies as in- 
cluding all galaxies listed in the Asiago Catalogue with 
Vr ^ 14,000km , not classified as Seyfert (including un- 
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Table 1. Distribution of of Supernovae 



SN Type 


Number of Supernova 


Events'' 




AU^ 


V,. < 14000 






Normal 


Seyfert 


Unclassified 


661 


36 


1 


I 


77 


3 


1 


la 


565 


192 


14 


Ib-Ic 


81 


48 


6 


II 


374 


210 


6 


St 


0.216 


0.229 


1.000 


Pec 


6 








Total 


1764 


489 


28 



'In the time lapse from Jan. 1, 1986 to Aug. 23, 2000 

''AH supernovas listed in the Asiago catalogue, since 1885 until 

August 23, 2000. 

" f= N(Ib/c)/N(II), as reported in the two preceding rows. 

certain Seyfert types, coded as "S"), Hll (code H2) or LIN- 
ERs (code S3) in the Veron-Cetty and Veron lists. The 
Seyfert sample comprises all Seyfert galaxies as above but 
coded as SI, S2 or Sla; in the catalogue (see Table 1 for 
sample sizes). 

We restrict our comparison to the rate of type Ib-Ic and 
II supernova, since the control ttme Tc is approximately the 
same for both type Ib/c and type II supernovce. The differ- 
ence in brightness and hght curve shape between type Ib-Ic 
an d type II supernovae wo uld translate in a /S.Tc/tc ~ 15 
% ( [Cappellaro et al. 1993 ), which is marginal with respect 
to the magnitude of the effects we are considering (a factor 
« 4; see §m). Even if Seyfert galaxies and normal galaxies 
have been monitored over different times, and with different 
techniques, the ratio N(Ib/c)/N(II) should a/most reproduce 
the true ratio of supernova rates for type Ib/c and type II 
(but clearly not each supernova rate in supernova units!). In 
other words, we have SNR(Ib/c)/SNR(II) = N(Ib/c)/(0.01 
Tibja) X 0.01 r7//N(II) « N(Ib/c)/N(II), where r is in cen- 
turies, SNR is expressed in number of supernova per year, 
and rc(II)~ rc(Ib/c). Selection effects, described in the next 
section, will operate in the same way in both galaxy samples. 

The Seyfert sample is dominated by Seyfert 2 galaxies, 
and at any rate by very low luminosity AGN. The average 
absolute magnitude of the Seyfert 2 sample is Mb ~ —19.1 
(taking the data of the Veron-Cetty and Veron Catalogue) 
The Six nuclei are, on average, of similar luminosity {Mb ~ 
— 19.3). In complete samples, Seyfert 1 nuclei are on average 
a factor 100 brighter than Seyfert 2 (see also §0). 



3 RESULTS 

3.1 The Distribution of Core Collapse Supernova 
Types 

Table 1 reports the number of supernova counted in each 
sample, with the restrictions on redshift and on discovery 
time motivated above. 

The ratio S = N(Ib-Ic)/N(II) for normal galaxies is 
RiO.23, very close to the values obtained from the ratio of ab- 
solute supernova rates in SNu: 0.18 for SOa-Sb g alaxies, and 
0.27 for Sbc-Sd galaxies ( Cappellaro et al. 1997 ). The ratio 
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Figure 1. Distribution of supernova distance normalized to the 
photometric diameter D25, for normal (upper panel) and Seyfert 
(lower panel) galaxies with at least one supernova listed in the 
Asiago supernovae catalogue. Supernova offsets and photometric 
diameters are as reported in the Asiago catalogue. Abscissa scale 
is ratio between deprojected supernova distance and photometric 
radius; ordinate scale is number of supernovas. Thick solid line: 
all supernovae; dotted lines: type Ib/c supernovae, dot-dashed line: 
type II supernova. 



S is much larger in the Seyfert than in the normal galaxies 
samples, by a factor « 4. Scc,Sy = N{Ib/Ic) /{Ncc,~ 0.5 S> 
Sr^orm ~ 0.19, whcrC Nee = N(Ib) + N(Ic) + N(II) is the 
number of core-collapse supernovae. Assumption of binomial 



statistics, where the probability is p=Sc 



-I, yields an ex- 



cess of type Ib-Ic supernova in the Seyfert sample significant 
to a confidence level of ^ 2a (N(Ib/c) expected in the Seyfert 
sample is 2.3±1.3). However, this simple estimate may not 
be fully satisfactory, since normal galaxies and Seyfert sam- 
ples may be influenced by factors affecting the probability of 
discovery as well as the intrinsic rate of supernova explosion 
(see 



3.2 The Radial Distribution of Supernovae 

Fig. |l] shows the distribution of the supernova de-projected 
distance d* normalized to the standard isophotal diameter 
D25, for normal and Seyfert galaxies, computed under the as- 
sumptions that all galaxies can be considered flat disks. Su- 
pernova in Seyfert galaxies show a tendency to occur closer 
to the nucleus. A Kolmogorov - Smirnov test suggests that 
the distributions are different to a statistical sign iflcance of 
0.999. This result has been known for some time (Petrosian 



h Turatto 199C; Kazarian 1997) from the analysis of super- 



nova data in different samples. 



3.3 Significance of Results: Is there any Real 
Excess of type Ib-Ic Supernovae in Seyfert 
galaxies? 

The ratio Snorm has been determined from 210 type II and 
48 type Ib-Ic supernovae; therefore it is possible to safely 
use the binomial statistics to write Snorm ~ 0.23 ± 0.05 
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Figure 2. Distribution of the ratio S= N{Ib/Ic)/N(II) for 20000 
bootstrap replications of the Seyfert sample by pseudo-control 
samples of the same size (see text). Thick solid lines: without re- 
distribution of unknown and type I supernovae; thin solid line: 
after redistribition of unknown and type I supernova. We con- 
sidered 25 galaxies with Vr < 14000km s~^ with supernovae dis- 
covered in or after 1985. The vertical dot-dashed line marks the 
value found for the Seyfert sample, as reported in Table 1. 
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(at 1 a confidence level). To properly estimate the statis- 
tical significance of the excess of type Ib-Ic supernova in 
Seyfert galaxies, we used a "bootstrap" method. We gener- 
ated a large number (~ 10*) of synthetic data sets matching 
the Seyfert galaxy sample. The synthetic data set were ex- 
tracted randomly from our normal galaxies samples, with 
the same number of galaxies as in the Seyfert samples, and 
with the same distribution of any parameter influencing ei- 
ther the probability of discovery of supernovae, or the intrin - 
sic rate of supernova production (Cappellaro et al. 1993): 



(1) galaxy inclination, (2) D25, (3) morphological type, (4) 
radial velocity. "Matching" samples were considered all syn- 
thetic samples for which the four distributions were all not 
statistically different from the distributions for the Seyfert 
sample at a 2a confidence level. In other words, since we 
could not eliminate selection effects from our Seyfert sample, 
we reproduced them in a similar fashion in a large number 
of "pseudo-control samples." 

In all cases save the morphological type, the statistical 
difference between the parameter distribution of the Seyfert 
galaxy sample and of same-size samples of normal galaxies 
was assessed by a Kolmogorov-Smirnov test. In the case of 
morphological type, a test was applied, after rebinning of 
T values into four bins for early type galaxies (E-SO), early 
spirals (SOa-Sb) late type galaxies (Sc and later), and un- 
known or peculiar. We computed the number of supernovae 
for each synthetic data set, and then we counted which frac- 
tion of synthetic data sets had supernova numbers exceeding 
the number from the real data. That fraction has been as- 



sumed to be our statistical significance ( Efron fc Tibshirani 
199CI ^ 

Fig. Q the thick solid line shows the distribution of S 
for the pseudo-control samples. The dot-dashed line marks 
the value found for the Seyfert sample. Only in ~ 1 % of 



Figure 3. Upper panel: expected ratio S of SNIbc to SNII rates 
as a function of the binary fraction. A Salpeter IMF and a con- 
stant star formation rate in the last 100 Myr have been adopted. 
Several cases are considered as detailed in the text. Middle panel: 
IMF effects. Ratio S in model c), for different IMF exponents x, 
as labelled in the figure. The Salpeter IMF corresponds to x=1.35. 
Lower panel: age effects. Ratio S in model c) for different star- 
burst ages (assuming solar composition), as labelled in Myr. 



the cases a larger value of the ratio N (Ib/c) /N (II) is found, 
confirming that a realistic estimate of the significance of the 
result is ^ 2a. The thin line marks the distribution of su- 
pernova after supernovae whose classification was unknown 
and or of generic type "I" h ave been re-assigned (cons is- 
tently with the approach of (Petrosian & Turatto 199C)). 
From Fig. ^ we infer that the redistribution for unclassi- 
fied supernova types has a negligible effect on the statistical 
results. 



4 INFERENCES ON STAR FORMATION 
4.1 Normal Galaxies 

The interpretation of the ratio 5 for the "average normal" 
galaxy is straightforward, if we consider most likely that Ib/c 
supernovae progenitors are Wolf Rayet stars, who se hydro- 
gen envelope has been fully lost (Filippenko 1998, and ref- 



erences therein). The latter stars may resul t from the evolu- 



tion of the most massive stars (M>4OM0)(Bressan, Chios! 



Bertelli 1981 ; Bressan 1994 ) or from less massive primary 
component of binary systems after the mass transfer (eg. 
Table 2 of (Wellstein & Langer 1999)). Assuming a constant 
star formation rate during the last 100 Myr, and adopting 
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the same initial mass function (IMF) for single stars and 
binary components, the ratio S is simply given by 



S = 



(1) 



where ih Iti the Itiitial triiirfs 6f th^i Sitar, th^; Itidex i ruiats t6 
all channels leading to type Ib/c supernova, j to the forma- 
tion of type II supernovae and W's are relative weights that 
depend on the overall binary fraction and on the relative 
frequency of case A, B and C in the binary evolution. 

The expected variation of the ratio <S is shown in the 
upper panel of Figj^, as a function of the binary fraction. 
To this purpose we have adopted a Salpeter IMF up to 
120 M0, and the values of m' and m-' were taken from Ta- 
ble 2 of (Wellstein fc Langer 199i^); the latter authors also 



assume that all stars more massive than 40 M0 explode as 
type Ib/c supernova, while single stars and secondary stars 
less massive than 40 M0 explode as type II supernova. We 
have considered several cases: (a) all binaries evolve as Case 
A (dotted line); (b) all binaries evolve as Case B (dashed 
line); (c) Case A and B share the same relative weight (solid 
line); (d) as in (c) but in the case that only binary mem- 
bers explode as type Ib/c supernova (lower dot long-dashed 
line); (e) as in (c) but in this case the initial mass for Ib/c 
supernova has been lowered from 40 M0 to 25 M0 (up- 
per dot short-dashed line). The horizontal line represents 
the observed average value for the normal galaxy sample 
(CS). We did not account for Case C evolution because it is 
considered to be a very rare event and as such it will only 
slightly modify the picture depicted in Fig.^. The fraction 
S does not appear much affected by the precise details of 
the binary evolution. For normal spirals, the observed ra- 
tio is consistent with a scenario where the progenitors of 
type Ib/c supernova are Wolf Rayet stars, either evolved 
from stars initially more massive than 40 Mq or from less 
massive primary components of binary systems. The data 
support a universal Salpeter IMF and indicate a binary frac- 
tion of about the 50%. A significantly higher binary fraction 
(about 75%) is required if one insists that only binary evo- 
lution give rise to type Ib/c supernova. As discussed by 
Wellstein & Langer (1999), uncertainties in the mass- loss 
rate, convection treatment and chemical composition may 
affect both the upper mass limit (40 M0) leading to type 
II supernova in single stars, and the mass transfer phase in 
binary evolution. Fortunately, uncertainties leading to dif- 
ferent type lb and type Ic supernova events are not relevant 
here, because we are considering their combined frequency. 



4.2 Seyfert Galaxies 

From Table 1 and Fig.^, we also infer that things could be 
marked ly different in Sevfert galaxies, where the observed 



of the radial distribution of Hll regions in normal, LINER 
and Seyfert galaxies agrees with a large number of detailed 
spectroscopic and imaging studies of Seyfert galaxies. Con- 
sistently, the radial distribution of supernova appears to 
be more strongly peaked in Seyfert galaxies than in normal 



galaxies in several samples (Kazarian 1997; Petrosian & Tu- 



ratto 199C) as well as in our own. However, enhancement 



of star formation rate in the inner galactic regions cannot 
explain the relative excess of type Ib/c supernova found in 
Seyfert galaxies. 

The observed radial gradie nts of chemical composition 
( |Vila-Costas fc Edmunds 1992| ) and the central concentra- 
tion of type Ib/c supernova in Seyfert galaxies (Fig. ^) im- 
ply that Ib/c supernova could be typically produced in 2 
Zq regions in Seyfert and in Zq regions in normal galax- 
ies. Because mass loss increases with metallicity and, con- 
sequently, the minimum initial mass for single Wolf Rayet 
production decreases, S could be a function of metallicity. 
As shown by model (d) in the upper panel of Fig.^, S could 
increase by the 100% if the minimum mass for Ib/c produc- 
tion decreases from 40 M0 to 25 Mq. However, this is not 
very likely because mass-loss is also sensitive to the lumi- 
nosity and the latter is a strong function of the initial mass 
( [Bressan 199^ ). 

Another possible explanation is that star formation in 
circum-nuclear and disk regions of Seyfert galaxies is a high 
binary fraction. If we assume that all young stars are in bi- 
naries and undergo mass transfer during their life, we would 
obtain 5~0.6. This value, though lower than the observed 
one, would double the expected relative frequency of binary 
stars with respect to the normal galaxies. Inspection of Table 
1 shows that the observed ratio N(Ia) /Ncc is 0.63 for the nor- 
mal galaxies and rises to 1 in the Seyfert sample. Would this 
ratio be free from selection effects, it would strongly indicate 
a larger probability of assembling binary systems. However, 
the magnitude at maximum and light-curve of type la super- 
nova is very different from that of core-collapse supernova. 
This makes the comparison of the ratio N(Ia) over number 
of core-collapse supernova meaningless if not accompanied 
by a careful examination of the relative control times which, 
unfortunately cannot be done here. 

An appealing alternative is constituted by a "top- 
loaded" IMF. The high efficiency at which gas is repro- 
cessed into stars is often taken as an indication in favor 
of a top-heavy initial m ass function in compact Starbursts 
e.g. ( Franceschini 1999). The middle panel of Fig. ^ de- 
picts the effect of changing the IMF exponent x, where 



-(l + x) 



and x=1.35 corresponds to the Salpeter 



IMF. The value of S increases as the IMF becomes flatter. 
The observed value in Seyfert galaxies would indicate a very 
flat IMF, above 8 AIq. We must however stress that conclu- 
sive evidence in favor of an IMF slope significantly flatter 
than the Salpeter one has not been found in exter nal galax- 
ies, not even in Wolf-Ra y et and Starburst galaxies ( schaerer 



ratio is S~l. Several Seyfert galaxies, of which a sizable 



fraction belongs to strongly interacting systems (?; Dultzin 



Contini, fc Kunth 1999; Leitherer(1999) ) 



Hacyan et al. 1999, ~ 2/3 of Seyfert 2), show significant, or 
even dominant (over the active nucleus luminosity) circum- 
nuclear star formation. A recent Hq narrow b and survey 
( Gonzalez-Delgado, Heckman, & Leitherer 2001) found that 
about 1/2 of Seyfert 1 and Seyfert 2 show maximum value of 
Hll region emission in the circum-nuclear regions. This study 



The necessity of a top heavy IMF disappears if, in 
Seyferts, we are observing a very young burst of star for- 
mation. To illustrate this effect we show in the lower panel 
of FigJ^ the run of the ratio 5 with the binary fraction, for 
different burst ages and a constant SFR. If the starburst is 
sufficiently old, let us say of age r* ^ 50 Myr, then there 
may exist evolved stars with any mass between the upper 
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IMF limit and the critical mass M~8 Mq for core collapse 
supernovae. But if the burst is younger, then all stars with 
initial mass below a given value - M(r) - are still on the 
Main Sequence and only stars more massive than M(r*) 
can explode. Different curves in the lower panel of Fig. ^ are 
labeled by the duration of the burst, in Myr. A Salpeter IMF 
is consistent with the observed Sseyf if the star formation 
episode is younger than « 13 Myr. This value depends only 
slightly on the assumed chemical composition. 



5 CONCLUSIONS 

We found that normal galaxies show Snorm ~0.27, a value 
consistent with the ratio of absolute supernova rates in SNu. 
With the caveat of our still-limited sample size, Seyfert 
galaxies show a peculiar distribution of supernova types, 
with higher frequency of type Ib/c supernovae {Sseyf ~ 1) 
than non-active galaxies. These finding are consistent with a 
"normal" supernova rate related to secular star formation as 
far as non-active galaxies are concerned. A scenario where 
Wolf Rayet stars are produced by the most massive stars 
(>40Mo) and by less massive (8M0<M<4OM0) primary 
stars of binary systems that soon or later undergo mass- 
transfer, reproduces fairly well the observed ratio, without 
particular assumptions on the IMF. 

A large type Ib/c rate with respect to type II super- 
novae, as found for Seyferts, is a new result. We consider 
several explanations as possible, namely a lowering of the 
upper mass limit for type II precursors due to a higher metal- 
licity, a very high binary fraction, a top heavy IMF and a 
young age of the burst of star formation. While metallicity 
is likely to act in the observed direction, because it enhances 
the mass-loss rate, the amplitude of the observed effect may 
not be explained. 

Even if the S value found for Seyfert galaxies needs 
confirmation with better statistics, it is interesting to note 
that an high S could be associated with the enhanced SFR 
in the circum- nuclear and regions of galaxies hosting Sey fert 
nuclei, especially of type 2 nuclei (Dultzin-Hacyan(1995)). A 
young age of the burst of star formation - probably triggered 



by interaction wit h a massive companion galaxy (Dultzin 
Hacyan et al. 199!; ) - is a likely cause of the high 5 observed 



among Seyfert galaxies, as it requires the minor number of 
assumptions. This interpretation is also consistent with an 
evolutionary sequence leading from IR-luminous Starbursts, 
to Seyfert 2 and eventually to Seyfert 1 galaxies, with the 
youngest Seyfert 1 being the ones of lowe st intrinsic power 
and/or the ones mos t heavily obscured ( |Gu, Maiolino, & 
Dultzin-Hacyan 200 1|). 
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